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Abstract 

We present a systematic approach to implementation of basic quantum logic gates operating on 
polar molecules in pendular states as qubits for a quantum computer. A static electric field prevents 
quenching of the dipole moments by rotation, thereby creating the pendular states; also, the field 
gradient enables distinguishing among qubit sites. Multi- Target Optimal Control Theory (MTOCT) 
is used as a means of optimizing the initial-to-target transition probability via a laser field. We give 
detailed calculations for the SrO molecule, a favorite candidate for proposed quantum computers. 
Our simulation results indicate that NOT, Hadamard and CNOT gates can be realized with high 
fidelity, as high as 0.985, for such pendular qubit states. 
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I. INTRODUCTION 



Quantum computers take direct advantage of superposition and entanglement to perform 
computations. Because quantum algorithms compute in ways which classical computers can- 
not, for certain problems they provide exponential speedups over their classical counterparts 







8|. That prospect has fostered a variety of proposals suggesting means to implement such 
a device js, [ol [lO|. DeMille has detailed a prototype design for quantum computation using 
ultracold polar molecules, trapped in a one-dimensional optical lattice, partially oriented 
in an external electric field, and coupled by the dipole-dipole interaction |^|. This offers a 
promising platform for quantum conaputing because scale-up appears feasible to obtain large 
networks of coupled qubits |ll|-ll3|. 



In previous work, we focused on entanglement and on consequences of using a strong 
external electric field with appreciable gradient, required to prevent quenching of the dipole 
moments by rotation and to enable addressing individual qubit sites The molecules 

were represented as identical, rigid dipoles a fixed distance apart and undergoing pendular 
oscillations imposed by the external electric field. We determined the dependence of the 
entanglement of the pendular qubit states, as measured by the concurrence function, on 
three unitless variables, all scaled by the rotational constant. The first specifies the Stark 
energy and intrinsic angular shape of the qubits; the second specifies the magnitude of the 
dipole-dipole coupling; the third variable specifies the thermal energy. Under conditions 
deemed amenable for proposed quantum computers, we found that both the concurrence 
and a key frequency shift, Au, that has a major role in logic gates, become very small for 
the ground eigenstate. In order that such weak entanglement can suffice for operation of 
logic gates, the resolution must be high enough to detect the Au shift unambiguously. 

For diatomic molecules, the Stark effect is second-order, therefore a sizable external elec- 



tric field is required to produce the requisite dipole moments in the laboratory 



subsequent study, we examined symmetric top molecules as candidate qubits 15||. Symmet 



rame. In a 



ric top molecules offer advantages resulting from a first-order Stark effect, which renders the 
effective dipole moments nearly independent of the field strength. That permits the use of 
a much lower external field strength in addressing sites. Moreover, for a particular choice of 
qubits, the electric dipole interactions become isomorphous with NMR systems. 

Here we study further aspects of how to implement a set of basic quantum gates for pen- 
dular qubit states of polar diatomic or linear molecules. We apply the Multi- Target Optical 
Control Theory (MTOCT) iJEl to design laser pulses that enable resolving and inducing 



2 



transitions between specified states of tlie qubit system. Tliis approacli lias been previously 
employed to study optimal control for elements of quantum computation in molecular sys- 



19 



27| . Our use of pendular 



tems, using as qubits vibrational or rotational states 
qubit states incorporates more fully effects of the external electric field and thereby simplifies 
the gate operations. 

Section II specifies the Hamiltonian defining the pendular qubits, as well as fundamental 
aspects of MTOCT. In Sec. Ill we present simulation results using MTOCT to obtain 
optimized laser pulses for realizing NOT, Hadamard and CNOT logic gates; those gates with 
the addition of the phase gate vr/S provide the basis for universal quantum computation [28]. 
Section IV discusses strategies to contend with cases in which the Au shift is zero or becomes 
too small to resolve. 



II. THEORY 

A. Eigenstates for Polar Molecules in Pendular States 

The Hamiltonian for an individual trapped polar diatomic or linear molecule in an external 
electric field e can, for our purposes, be reduced just to the rotational kinetic energy and 
Stark interaction terms [14]: 

rCs = B ■ - fiecosO (1) 

This represents a spherical pendulum: BJ^ is the rotational energy, with B the rotational 
constant, and fi the permanent dipole moment; 6 is the polar angle between the molecular 
axis and the external field direction. At the ultracold temperatures that we consider, the 
translational kinetic energy of the trapped molecules is very small and nearly harmonic within 
the trapping well, so the trapping energy is nearly constant and hence is omitted. Interactions 
involving open shell electronic structure or nuclear spins or quadrupole moments are also 



omitted (but could be incorporated in familiar ways 29l|). The eigenstates of resulting 
from mixing of the field-free rotational states by the Stark interaction, are designated as 
pendular states. As proposed by DeMille jg], the qubits |0) and |1) are chosen as the two 
lowest M = pendular states, with J = and 1, respectively. Here the J notation (tilde- 
hat) indicates it is no longer a good quantum number, due to the Stark mixing. However, 
M remains good as long as azimuthal symmetry about e is maintained. The qubits thus are 
superpositions of spherical harmonics, 
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|0) = 5^a,-y,,o(^, y.); \l) = J2bry^,o{9,^) (2) 
j j 
Adding a second molecule into the trap, identical to the first but at distance ri2 from 

it, introduces in addition to its pendular term the dipole-dipole interaction, Vdd- Averaging 

over the azimuthal angles [14], which for M = states are uniformly distributed, reduces 

the dipole-dipole interaction to 

Vdd = ■ COS 01 ■ COS 02 (3) 

where fi^ = (1 — 3 cos^ a), with Q = ^Vrfj. As depicted in Fig. [T], a is the angle between 
the array axis and the electric field direction e; 0i and 02 are the polar angles between the 
dipoles and the field direction. 

To exemplify logic gate operations, it is sufficient to consider just two molecules. In the 
basis set of qubit pendular states {|00) , |01) , |10) , |11)}, the two-molecule Hamiltonian can 
be expressed as Ktot = 'Hsi + "^52 + K/^, with 



U 



SI 



dd 




(g) I2; = I2 ® 



/^i (-11 



(4) 



Here Wq and Wi are eigenenergies of the pendular qubits for states |0) and I2 is the 
2x2 identity matrix. Co, Ci are the expectation values of cos^ in the basis of |0) and |1); 
while Cx indicates the transition dipole moment between |0) and The matrix elements 
thus are defined as: 



Co = (O|cos^|0); Ci = (l|cos^|l); C^ = (0 |cos ^| 1) (5) 

In Eq. [3l primes are used to indicate that the external field strength differs at the location 
of the two molecules, as required to distinguish between the qubit sites. The entanglement 
of the pendular qubit states of Hfot was evaluated in ref jl^. From numerical results, simple 
approximate formulas were obtained that provide the concurrence and the key frequency 
shift, Acu, in terms of two unitless reduced variables: x = ^^^/B and ^q^/b. When ^^q/b ^ 1, 
the usual case, the concurrence is proportional to Aw, which is given by 

Aa; = |fi„|(Ci-Co)(C;-C;) (6) 
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Figure [2] plots this relation, wherein ^'^/n^ depends only on x and x' — x. 



B. Multi- Target Optical Control Theory (MTOCT) 



Much attention has recently been devoted to appl ying Op timal Control Theory for el- 
ements of quantum computation in molecular systems 16l-|l9|. The basic idea is to design 
laser pulses which allow manipulation of transitions within each qubit separately. For im- 
plementing basic quantum gates, the aim is to achieve large transition probabilities with 
the correct phase from a specific initial state into a final target state by application of an 
external laser field while minimizing the laser energy. For our case, we can construct the 
following MTOCT objective function, 53, which needs to be maximized 23] : 
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X J l^fk it) 




i^^k {t) )dt)) (7) 



where z is the total number of targets and for qubits is given by z = 2^ + 1, where 2^ 
is the number of input-output transitions in the gate transformation and the supplementary 
equation is the phase constraint. Thus, for the two dipole system, z = 'b. Here, tpi^ is 
the wave function of the /c-th target, driven by the laser field E (t), with initial condition 
V'ifc (0) = (pik- Whereas ipf^ is the wave function of the A;-th target driven by the same laser 
with final condition tpfk (T) = (pfk- Thus, the first term on the right-hand side represents the 
overlap between the laser driven wavefunctions and the desired target states. In the second 
term, E (t) is the laser intensity, S (t) = siri^ {^^/t) is the laser envelop function, which 
guarantees the experimentally appropriate slow turn-on and turn-off laser pulse envelope 
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30] . T is the total duration time of the laser, is a positive penalty factor chosen 
to weight the importance of the laser fiuence. The last term denotes the time-dependent 
Schrodinger equations for wave functions ipik{t) and ipfk{t) with l-Ltot the Hamiltonian of the 
system. 

Requiring (5^3 = 0, specifies the equations satisfied by the wave function, the Lagrange 
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multiplier, and the optimized laser field 23| : 

d 



th—J^k (t) = {Htot -fi-E{t)} (t) 
at 

d 

^ifc (0) = (pik, i'fk (T) = 4>fk 

and k = 1 ■ ■ ■ z 



h ■ ao 



k=l 



{ipfk (t) \cos9i + cos92\i)ik {t))} 



(9) 



In order to examine the performance of the optimized laser pulse, we evaluate two factors 
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31 



32|: the average transition probability, given by. 



(10) 



k=l 



and the fidelity, given by 



Y.{i'^k (T)|$/,) 



k=l 



The average transition probability involves only the overlap of ipi]^ (T) and but does 
not reflect the difference of phase information between the laser driven final state ipik (T) and 
the designed target state ^jk- Since the phase information is very important for quantum 
logical gates, we use only the fidelity parameter to assess our simulation results. 

III. SIMULATION RESULTS FOR POLAR DIATOMIC MOLECULES 



A number of diatomic polar molecules offer properties suitable for a quantum computer 

m 
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23 



25|. For our numerical study, we chose SrO, for which the dipole moment /i 



8.9 Debye and rotational constant B = 0.33 cm~^ [33]. Since we consider trap temperatures 
in the microkelvin range, with ^BTj^ 10~^, thermal excitations are negligible. To specify 
the pendular states and other properties requires assigning the external field strengths at 
the sites of the two molecules and the distance between them. We used field strengths (in 
scaled units) of /ie = 2 and 3, corresponding to e = 4.4 and 6.6 kV/cm, respectively, at 
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the two sites. Initial 
an optical lattice 



we took ri2 = 500 nm, a typical spacing for molecules trapped in 
. We set the angle a = 90° (cf. Fig. [1]), the usual experimental 
choice. Then ^^/b = 9.7 x 10~^. The corresponding pendular eigenstate reduced energies are 
e,/b = -1.65, 1.19, 1.92, 4.77 and the cosine matrix elements are Co = 0.480, Ci = -0.208; 



C'q = 0.579 and C'l = —0.164. Thus, from Eq. [6] we obtain Au = 51 kHz. As seen in Fig. 
|2l Au varies only modestly with the field strength in the range i^^/b = 2 — 5, considered 
optimum , but Au is directly proportional to the dipole-dipole coupling strength, fi^. 

At present, it remains an open question whether, in the presence of line broadening 
induced by the static external electric field, adequate resolution can be obtained to resolve 
unambiguously a frequency shift of only 50 kHz 1^. Such a small Atu is also a severe 
handicap for our theoretical simulation of laser-driven logic gates. For instance, the laser 
pulse duration 26| required for realizing the CNOT gate is r = io^aw. Hence a frequency 
shift so small as 50 kHz requires that the laser pulse duration is at least 31/iS. Recently, Zaari 
and Brown studied the effect of laser pulse shaping parameters about the fidelity in realizing 
the quantum gates. They proved that the amplitude variation and frequency resolution plays 
the important role in the fidelity 3J|. We will explore the impacts of those two coefficients 
in the future work. That is much longer than the self-evolution period of the system, about 
29 ps. If we set the computation time step at 0.25 ps, a single simulation run would need 
1.25 X 10^ steps. Our MTOCT calculations involve many iterative runs; e. g., for a CNOT 
gate about 440 iterations. With current computers, such a calculation would be daunting: it 
would need about 700 GB of RAM storage to perform and each iteration would take about 
2.8 days, in total nearly 3.5 years for standard 17 core. 

To make the calculation feasible, we reduced the spacing between the dipoles ten-fold, 
which increases Q by 1000-fold, and thus Au = 50 MHz and the laser pulse duration shortens 
to 33 ns. Then simulation runs have ~ 10^ steps, the RAM needed shrinks to 700 MB, and 
the computation time to about 4 minutes per iteration, so ~ 30 hours for 440 iterations. 
The reduced dipole-dipole spacing, which becomes only 50 nm, actually corresponds to the 
range recently proposed for plasma-enhanced, electric/electrooptical traps, for which the trap 
frequencies can exceed 100 MHz 35M37|. and might be attainable in an optical ferris wheel 
device 38 



Reducing the spacing so markedly is not considered practical, however, because it 
would strongly foster inelastic, spontaneous Raman scattering of lattice photons and hence 

BE 



induce unacceptably large decoherence (5|, l6|, [39 1. The resort to taking ri2 unrealistically 
small was done reluctantly but enables us to illustrate the general utility of MTOCT applied 
to quantum logic gates. 
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In the simulations, the time evolution of (t) and \1/ fk (t) is calculated from Eq. [8] by 

lod, using time steps of 0.25 ps. The penalty factor is 
23] ■ The penalty factor is used to minimize the fluence of 



the fourth-order Runge-Kutta met 
set as 5 X 10^, the same as in Ref. 
the external fields 



25| . For the optimized laser pulse E(t) from Eq. [HI we adopted a rapidly 



convergent iteration using a first-order split-operator approach 40|. The maximum iteration 
number was set at 600. In addition to the usual four-qubit basis set, {|00) , |01) , |10) , |11)}, 
we also included the phase correction, introduced into the MTOCT approach by Tesch and 
de Vivie-Riedle [17]. In ref. [23j], Mishima and Yamashita pointed out that the purpose of the 
phase constraint is preventing each state evolving to different phases, which can provide the 
correct quantum logical gates. Recently, Zaari and Brown pointed out that align the phase 
of all qubits appropriately can lead to effective subsequent quantum gates (laser pulses) l^^. 

In Table m we show both the initial and target states for the NOT gate for the two dipoles. 
The optimized laser pulse for NOT gates applied to dipole 1 and dipole 2 separately is shown 
in the upper panel of Fig. [31 The pulse for the NOT gate applied on dipole 1 was obtained 
after 146 iterations with the converged fidelity difference of 1.26 x 10^^; the pulse for dipole 2 
took 76 iterations with the converged fidelity of 0.985. In our simulation, the total iteration 
number is decided by two main effects, the fidelity and the difference between the fidelity 
and the fidelity of the previous step. If the fidelity is greater than 0.9 and the difference 
is smaller than 10~^, the iteration will stop. Both pulses have similar maximum intensity, 
about 1.2 kV/cm. In order to verify the performance of this pulse as well as the time evolu- 
tion of the system within the laser pulse, we chose an initial state cos (|) |01) + sin (|) |10) 
and examined the population evolution. For an ideal pulse, the final state for the action 
of a NOT gate on dipole 1 should yield a final state of cos (|) |11) + sin (|) |00). Corre- 
spondingly, the final state for a NOT gate on dipole 2 should be cos (|) |00) + sin (|) |11). 
The population evolution due to this pulse assistance is shown in the lower panel of Fig. |3l 
Before achieving the final state, the population of each state oscillates a number of cycles. 
The NOT pulse for dipole 1 produced the converged population for state |00) to be 0.756 
and that of |11) to be 0.240, while the pulse for dipole 2 yielded populations of 0.263 and 
0.723 for states |00) and |11), respectively. As this final yield result approaches the ideal 
case (cos^60° = 0.75; sin^ 60° = 0.25), we conclude that the optimal laser pulse drives the 
system from an initial state to a target state according to Table [T] for NOT gates. 

The converged laser pulses by which to realize the Hadamard gate for both dipoles are 
shown in Fig. H] and the initial and target states are given in Table [Tll The fidelity for 
each pulse on each site is 0.944 and 0.902, respectively. The maximum intensity for the 
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pulse is around 1.8 kV/cm, slightly larger than that for NOT gates. We selected |00) as 
our initial state and plotted the time evolution in Fig. HI The population evolution of the 
Hadamard gate for dipole 1 is similar to the situation of the NOT gate (Fig. [3]) . The 
population of different states oscillates and finally yields 0.491 for |01) and 0.508 for |11). 
For the Hadamard gate on dipole 2, the population for |10) and |11) is almost zero during 
the entire evolution. But the population switches between |00) and |01) . This population 
oscillation commences at 6 ns and continues until the end of the pulse with a mean value 
of 0.5 and amplitude around 0.5. The final converged population is 0.529 for |00) and 0.470 
for |01). The population transfer, shown in Table IIT| going from an initial state to a target 
state, indicates that the design pulse successfully implemented the Hadamard gate. There 
is a small variation among all fidelity values we obtained in the simulation, which range 
from 0.902 to 0.985. This is due to the short duration of the laser pulse we are using in the 
simulation. The fidelity could be improved further by extending the duration of the laser 
pulse. 

Figj5] shows the converged laser pulse which performs the CNOT gate; the pulse is highly 
oscillatory with total duration time of 33 ns. The rapid oscillation is due to the short 
system evolution period. The maximum amplitude of the pulse is around 1.5 kV/cm, in 
the range easily achieved experimentally. The laser pulse was obtained after 441 iterations 
and yields fidelity of 0.975. If we use the same initial condition as that for NOT gates, the 
ideal final state should be cos (|) |01) + sin (l) |11). The simulated population evolution 
due to the laser pulse evolved is plotted in the lower panel of FigJS] The final population is 
0.0005, 0.2290, 0.0029 and 0.7675 for |00) , |01) , |10) and |11), respectively. The population 
evolution in Fig. [5] follows the initial and target results of Table IIIIl thus confirming the 
correct operation of the CNOT gate. 

In order to test, at least modestly, how the MTOCT approach responds to a change in 
pulse duration for the dipole-dipole system, we increased the spacing between the two dipoles 
to 75 nm. The frequency shift then becomes Aw = 14.6 MHz and the duration of the laser 
pulse is 110 ns. We carried out the simulations just for the CNOT gate. The optimized 
laser pulse is shown in Fig. O After 500 iterations, the converged fidelity is 0.90, which is 
slightly smaller than for the 50 nm spacing. Again, we tested one sample initial state. The 
population evolution is plotted in Fig. [61 As before, the population oscillation continued 
during the whole process and the final population obtained confirms the CNOT operation. 
This serves to indicate that the MTOCT approach is stable and provides a useful general 
means to implement logic gates for dipole-dipole systems. 
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IV. DISCUSSION 



We have applied the MTOCT methodology to pendular states for a pair of polar 
molecules (SrO) to determine the optimum laser pulse for implementing the NOT, CNOT 
and Hadamard quantum logic gates. Our results confirm that, for the conditions adopted 
(ri2 = 50 nm, resulting in Au ~ 50 MHz), a single laser pulse (with minimum duration 
~ 33 ns and amplitude < 2 kV/cm) suffices to operate these gates with high fidelity. How- 
ever, computational limitations (storage and time) did not permit us to treat conditions 
(ri2 = 500 nm, with ~ 50 kHz) that are considered congenial for experimental implemen- 
tation. 

This shortcoming is also manifest, in different ways, in two previous applications of 
MTOCT to assess laser-operated logic gates for polar molecules. Table |IV] compares our 
conditions with those studies. Both nominally emulated the design by DeMille with two 
ultracold trapped diatomic molecules entangled via dipole-dipole interaction. The version 
most akin to ours, presented by Bomble, et al 26|, considered a pair of NaCs molecules, with 
the external static field aligned along the intramolecular axis, ri2. Then a = 0° (rather than 
90° as in our case) and hence Qa = —2u). The qubits were taken as rotational states mixed 
by the Stark effect to second order (thus a fairly good approximation to the pendular eigen- 
states we used). The conditions adopted (ri2 = 300 nm, resulting in Au ~ 120kHz ) are 
considered suitable for experimental implementation. However, in carrying out the MTOCT 
computations for our system, a very large step size of 10 ps was used. That avoided entirely 
the storage and time limitations we encountered (with step size 0.25 ps). However, we found 
that replicate calculations using a 10 ps step size gave markedly irreproducible results for 
the optimal laser properties and population evolution of the logic gates under our simulation 
system. 



The other previous version, presented by Mishima and Yamashita 23|, |25| omits alto- 
gether an external static electric field, and takes as qubits the lowest "pure" rotational 
states (J = 0, M = and J = 1, M = 0) and lowest vibrational states of each molecule. 
Specific alignments of the molecular axis with respect to a laboratory fixed z-axis are con- 
sidered, but that is unrealistic because without an external electric field the molecular axis 
distribution is isotropic and the laboratory projections of the dipole moments vanish, as 
emphasized elsewhere . Also unrealistic from an experimental perspective is the choice 

of an extremely small distance between the molecules (ri2 = 5nm). That produces large 
entanglement but would induce severe decoherence {5], Q. The MTOCT treatment of logic 
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gates is nonetheless of interest, since the choice of "pure" rotational states as qubits results 
in Au = 0. Then, if the molecules are identical, transitions involving qubits on different 
sites cannot be resolved. In order to target individually the two sites, two different laser 
fields were used. In the MTOCT analysis, the separate laser pulses for gate operations then 
were much shorter and simpler than in our application using the same laser for both sites. 
For instance, to perform a CNOT gate using "pure" rotational states with two lasers took 



only 10'^ ps 23|, |25|, whereas our use of pendular qubits with one laser required up to 10^ 



ps (see Figs. OandlH])- Although the two-laser mode is theoretically inviting, it requires 
spatial resolution adequate for each laser to drive only one of the two molecules. That is not 
feasible unless the distance between the molecules is much larger than 5 nm. Another means 
to contend with Au = 0, or when it too small to resolve, has been exemplified in designs 
employing superconducting flux qubits 42|]. Again, that method requires spatial resolution 
sufficient to enable qubits on different sites to be driven individually. 

Taken together, the three studies of Table [IV] illustrate both the utility of MTOCT and 
the limitations imposed by present-day computational capability. Aptly, those limitations 
foster yearning for the arrival of a quantum computer. 
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Table I: The initial and target states of the NOT gate. The 
final fidelities for both NOT gates are 0.967 and 0.985, re- 
spectively. 





Initial State 


Target State 


i 


NOT gate for Molecule 1 


NOT gate for Molecule 2 


1 


|00) 


|10) 


|01) 


2 


|01) 


111) 


|00) 


3 


|10) 


|00) 


111) 


4 


111) 


|01) 


|10) 


5 


i(|00) + |01) + |10) + |ll)) 


i(E^<^-(|10) + |ll) + |00) + |01)) 


ie^'^-(|01) + |00) + |ll) + |10)) 



Table II: The initial and target states of the Hadmard gate. 
The yield fidelities are 0.944 and 0.902 





Initial State 


Target State 


i 


Hadamard gate for Molecule 1 


Hadamard gate for Molecule 2 


1 


|00) 




75(100) + |01)) 


2 


|01) 


;^(|01) + |ll)) 




3 


|10) 




^,m + m 


4 


111) 






5 


i(|00) + |01) + |10) + |ll)) 


^<E^<^-(|00) + |01)) 


^e^'^-(|00) + |10)) 



Table III: The initial and target states of the CNOT gate. 
The converged fidelity is 0.975. 



i 


Initial State 


Target State 


1 


|00) 


|00) 


2 


|01) 


|01) 


3 


|10) 


111) 


4 


111) 


|10) 


5 


i(|00) + |01) + |10) + |ll)) 


ie*'^-(|00) + |01) + |ll) + |10)) 
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Table IV: Comparison with other polar molecular systems 





Molecule ^ (D) R 


(nm) 


B [cmr^] 




X 


x' 


Atj (MHz) 


T (ns) 


Ref. [26] 


NaCs 


4.6 


300 


0.059 


6.7 X 10"^ 


2.0 


2.6 


0.12 


1.2 X 10^ 


This paper 


SrO 


8.9 


50 


0.33 


9.6 X 10^3 


2 


3 


51 


33 


Ref. [23, 25| 


NaCl 


8.4 


5 


0.22 


13 









2.6 


Ref. [23| 


NaBr 


8.2 


5 


0.15 


18 












15 







y 




/ 









At2 



Figure 1. The configuration of two polar molecules, /xi and /12 are permanent dipole moments for 
molecule 1 and 2; r\2 is the distance vector from molecule 1 to 2. e is the external electric field and 
a is the angle between ri2 and the external field. 9i (^2) is the angle between /ii {^l2) and e. 




0.1: : 

0.0: : 

0.0 0.2 0.4 0.6 0.8 1.0 

x'-x 

Figure 2. Ratio of frequency shift, Ao;, to the dipole-dipole interaction parameter, = 

(l — 3cos^a) /x^/r^2i as a function of reduced field strength, x = i^^/b and x' = p-^'/b at sites 
of the two dipoles. 
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NOT Gate for Dipole 1 NOT Gate for Dipole 2 



T ' 1 ' 1 1 I ' 1 ' 1 ' r 




Time (ps) 

Figure 3. Converged laser pulses for NOT gate. The upper panel is the laser pulses for realizing 
the NOT gate for dipole 1, on the left and dipole 2 on the right. The initial and target states are 
listed in Table HI The lower panel shows the evolution of all populations driven by NOT pulses. The 
left panel exhibits the population evolution under the NOT pulse for dipole 1. The initial state 
is cos (^) |01) + sin (^) jlO) and the final populations are 0.733 for |00) and 0.265 for |11). The 
right panel shows the population evolution via the NOT pulse for dipole 2. It has the same initial 
condition as the left panel and the converged population is 0.216 for |00) and 0.781 for [11). Both 
pulses perform the corresponding NOT gates nicely. 
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Hadamard Gate for Dipole 1 Hadamard Gate for Dipole 2 



T ' 1 ' 1 1 I ' 1 ' 1 ' r 




Time (ps) 



Figure 4. Optimized laser pulses for Hadamard gates. The left panel pertains to performing the 
Hadamard gate on dipole 1 and the right one to performing that on dipole 2. The lower panels show 
the population evolution under the Hadamard laser pulse. In order to exhibit the curves clearly, we 
chose |01) as the initial state; the final populations of |01) and |11) are 0.429 and 0.566, respectively. 
For the case of the Hadamard gate on dipole 2, the |00) qubit was chosen as the initial state. The 
populations at the end of the evolution are 0.528 for |00) and 0.470 for |01). The functions of these 
gates are realized very well by these two pulses. 
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CNOT Gate 




Time (ps) 



Figure 5. Optimized laser pulse for realizing the CNOT gate. The initial and target states are listed 
in Table IIIII In the lower panel, the population evolution is driven by the CNOT pulse. The initial 
state is cos (^) |01) + sin |10). After population oscillations due to the effect of the laser pulse, 
the final qubit populations are 0.00363 (|00)), 0.25807 ([01)), 0.00273 (|10)) and 0.73546 (|11)). 
The populations of state |10) and |11) are switched, which confirms the correctness of the converged 
laser pulse. 
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CNOT Gate 




Time (ps) 



Figure 6. Optimized laser pulse for realizing the CNOT gate when the distance between two dipoles 
is 75 nm. The initial and target states are listed in Table IIIII The optimized laser pulse, which is 
shown in the upper panel, lasts 110 ns. The lower panel shows the population evolution driven by 
the pulse. The initial state is the same as Fig. [S] The final converged populations for state |01) 
and 1 10) are 0.22 and 0.73, respectively. 
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